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. Abstract 

CN . The present work is an endeavour to connect the properties of tiny nearly massless 

objects with those of some of the most massive ones, the compact stars. 



O ! Since 1996 there is major influx of X-ray and 7 ray data from binary stars, one or 



both of which are compact objects that are difficult to explain as neutron stars since they 
contain a mass M in too small a radius R . The suggestion has been put forward that 
Q_|. these are strange quark stars (SS) explainable in a simple model with chiral symmetry 
restoration (CSR) for the quarks and the M, R and other properties like QPOs (quasi 
^ I periodic oscillations) in their X-ray power spectrum. 

^ ■ It would be nice if this astrophysical data could shed some light on fundamental 

properties of quarks obeying QCD. One can relate the strong coupling constant of QCD, 
to the quark mass through the Dyson-Schwinger gap equation using the real time 
rS I formalism of Dolan and Jackiw. This enables us to obtain the density dependence of 
I from the simple CSR referred to above. This way fundamental physics, difficult to extract 
from other models like for example lattice QCD, can be constrained from present - day 
compact star data and may be put back to modelling the dense quark phase of early 
universe. 
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A calculation was done for cold stars in |I| and important conclusions drawn from 
there about cliiral symmetry restoration in QCD when the EOS was used to get SS fitting 
definite mass-radius (M-R) relations [0, ^ ^. The empirical M-R relations were derived 
from astrophysical observations like luminosity variation and some properties of QPO-s in 
the X-ray power spectrum of these compact stars. The calculations are compared to these 
stars which emits the X-rays, generated from accretion from their binary partner. The 
different QPO-s show a correlation for not only these stars but black holes also so that 
one cannot explain these as being due to magnetic fields or other properties of the stars 
but some scale-dependent phenomenon in the orbits of the particles which are accreted 
p. One of the compact objects, the SAX J1808.4 - 3658 with period 2.49 millisecond, 
has been called the holy grail of X-ray astronomy Its discovery was anticipated for 
nearly 20 years because magnetospheric disk accretion theory as well as evolutionary ideas 
concerning the genesis of millisecond radio pulsars strongly suggested that such rapid spin 
frequencies must occur in accreting magnetic field neutron stars. 

The interesting point made by Dey et al. |1| is that starting from an empirical form 
for the density dependent masses of the up (u), down (d) and strange (s) quarks (q in 
short) given below, one can constrain the parameter of the form of this mass from recent 
astronomical dataQ: 

Pb 

Mi = rrii + Mq sech{iJ — ), i = u,d,s. (1) 

Po 

For three values of u = 0.286, 0.333 and 0.400 the q-star can have a sequence of 
masses obtained from the standard TOV equations for different choices of the central 
density. This is shown in Fig. (|l|). 

In other words the masses of stars in units of solar mass, {M/Mq), found as a function 
of the star radius R, calculated using the above eqn.(|l]), - produces constraints which 
enable us to restrict the parameter u. At high pB the q- mass Mj falls from its constituent 
value Mq to its current one m,. The other parameter Mq was taken to be 0.31 GeV to 
match up with constituent quark masses assuming the known fact that the hadrons have 
very little potential energy. The results are not very sensitive in so far as changing Mq 
to 0.32 GeV changes the maximum mass of the star from 1.43735 Mq to 1.43738 Mq and 
the corresponding radius changes from 7.0553 kms to 7.0558 kms. 

It is interesting to plot the up (u), down (d) and strange (s) quark masses at various 
radii in a star. This is done for a particular value of the parameter u = 1/3 already 

^ Pb = {pu + + Ps)l'i is the baryon number density, po — 0.17 fm~^ is the normal nuclear matter 
density, and ^' is a numerical parameter. The current quark masses rrii used in the following are 4, 7 and 
150 MeV for u, d and s respectively. 
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Figure 1: Variations of the mass with central densities pc (in units of 10^^ g/cm^), for 
three different values of the parameter v (= 0.286, 0.333 and 0.400) as found in [|^. 

discussed in [|T], ^ which gives a sequence of stars falling right into the allowed region 
of the M-R curve. Fig shows that the quarks do not have the constant current masses 
assumed in the bag model nor do they have the constituent masses of zero density hadrons. 
Upto a radius about 2 kms the quarks have their chiral mass but in the major portion of 
the star their masses are substantially higher. At the surface the strange q- mass is about 
0.278 GeV and the u,d q- masses ~ 0.13 GeV . 

Given the form for the mass eq. (|l]) we have to look for a formalism to calculate this 
mass in QCD. This can be done very conveniently using the real time formalism of Dolan 
and Jackiw ||5| (DJ in short) since here one does not encounter the problem of imaginary 
chemical potential. As is well known usually people work in the Euclidean space for 
lattice and other gauge invariant theories. But this excludes finite density calculations 
since that involves imaginary chemical potentials that makes the action unbounded from 
below. Using D J, the price one has to pay is that this is a gauge dependent formalism. We 
work in the well established Landau gauge used by . We believe in the physical nature 
of our results and wish to impress upon the reader the robustness of this physicality. In 
other words we argue that the modelling that we have to invoke does not take away the 
basic nature of the results that we obtain, namely that the strong coupling constant as 
decreases with increasing density and in principle this can be constrained by compact star 
data if one believes them to be SS. 
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Figure 2: Quark masses at different radii in the star. Notice that chiral symmetry is 
restored inside a radius of 2 kms but outside of this the masses are substantial. 



The Fermion propagator in the DJ formahsm is as follows : 

+ m) 2it6{p^ — m^) 



Sf{p) 



— rn? + ie exp[po — ^f/T] + 1 



(2) 



where T is the temperature and ep-, the chemical potential. There is a similar propagator 
for the gluon involving the Bose function instead of Fermi but with no chemical potential. 
One has to calculate the self energy of the quark with this propagator including a gluon 
loop of four momentum k in Fig(p|). 




Figure 3: The self energy of quarks with a gluon loop 
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Using 'y^D^i,'j'^ = — Si/k"^ and a colour factor 4/3 for the colour group SU(3) one 
getsi 

where we have neglected a Boson finite T-term which deos not contribute and / is the 
Fermi function ^ 

^^^^ " l + exp[iE-EF)/Ty 

In the rest frame, p = 0. The Fermi function becomes a step function for our case, since 
the temperature T = 0. On closing the contours in the lower half ko plane, first term 
reduces to 

d^k , 1 1 , Aq^m f dk 



(27r)3^2E[(m-E)2-fc2] 4^^^ gn^ J E 
while the second term becomes, on using the variable p — k —>■ q, 

.2 f ^(g^ - m(p)2) ^ -4c/2m r dq 2 

^'"^J {2ny{p-qy.{exp[{qo-fi)P] + l Svr^ i E exp[(go - + 1 

The first integral is logarithmically divergent and therefore must be set to renormalize 
the p = 0, T = 1/beta = quark mass. We do not use a high value of 01/^11 = = 0.75 
like since it gives a cut-off which is too small, of magnitude 1.29 GeV only. We prefer 
a lower = 0.55 and this also gives us a high cut-off ko = 2.86 GeV. We use 
m{p = 0) = 0.33 Gel/. 

Now we introduce the second term for finite p but T = so that the Fermi function 
reduces to a step function. 

as f^p dE , , 

m p = m - ^ / , 7 

Given our form for mlp), which is Mq of eq.([l|), dropping the current quark mass 
rrii, as is now evaluated for all densities. For small density its value increases but this 
is unnecessary for our model. We get very reasonable values of for p = 5po upwards 
as can be seen in Table 1 and also in Fig(^). The calculation serves double purposes, it 
shows that for the concerned densities our choice of the mass function, namely eq.(p, is 
reasonable and also as already stressed it gives a novel shape for the variation of as with 
density not easy to obtain otherwise. 



5 



Pb/po 




0.333 


0.286 


0.400 
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0.5093 


0.522 


0.5156 
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0.5077 


0.5025 


0.478 
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0.4919 


0.4739 


0.4384 
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0.4693 


0.4427 


0.4012 
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0.444 


0.412 


0.3676 
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0.4185 


0.3833 


0.3379 
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0.3939 


0.3571 


0.3119 
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0.3707 


0.3334 


0.2891 
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0.3493 


0.3121 


0.2691 
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0.3297 


0.2929 


0.2516 
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0.3117 


0.2758 


0.2361 
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0.2953 


0.2604 


0.2223 


17 




0.2803 


0.2465 


0.21 


18 




0.2666 


0.2339 


0.199 



Table 1: Variation of with density for three values of u 

The variation with u is very small at the surface but deep inside the star when the 
density increases the values of goes down substantially and the change assumes sig- 
nificance. The physicality of the results will enable us to extend our calculation now to 
T and apply to early universe where temperature of order 100 MeV is expected. 

In summary we have shown both in principle and in practice that the strong coupling 
constant can be derived at various densities once one knows the behaviour of the chiral 
symmetry restoration for the quark mass constrained from stellar data. The practical 
model that we have chosen uses an empirical form for chiral restoration at high density 
proposed and tested against various star properties through the allowed mass radius 
regions of the compact objects[|l], If. We can therefore urge the astrophysics data- 
analysts to pinpoint the mass and radius of compact objects like SAX J1808.8 to a greater 
precision to help efforts like ours. 



It is a pleasure to thank Dr. Arun Thampan for helpful discussions. The stimulation 
for the work came from a discussion with Prof. Donald Lynden-Bell. 
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Figure 4: The figure shows the smooth decrease in the value of the strong couphng 
constant with the increase of baryonic density. We have plotted three curves corresponding 
to the three values of the parameter z/ ( = 0.286, 0.333 and 0.040) of eq. (|l|). 
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